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Abstract — Though network coding is traditionally performed 
over finite fields, recent work on nested-lattice-based network 
coding suggests that, by allowing network coding over finite chain 
rings, more efficient physical-layer network coding schemes can 
be constructed. This paper considers the problem of communi- 
cation over the flnite-chain-ring matrix channel Y = AX + BZ, 
where X is the channel input, Y is the channel output, Z is 
random noise, and A and B are random transfer matrices. Tight 
capacity results are obtained and simple polynomial-complexity 
capacity-achieving coding schemes are provided under certain 
distributions of A, B, and Z, extending the work of Silva, 
Kschischang and Kotter (2010), who handled the case of finite 
fields. This extension is based on several new results, which may 
be of independent interest, that generalize concepts and methods 
from matrices over finite fields to matrices over finite chain rings. 

Index Terms — Lattice network coding, finite chain rings, ma- 
trix normal form, matrix channels, channel capacity. 



I. Introduction 

A MATRIX channel provides a useful abstraction for 
studying error control for linear network coding schemes. 
Transmitted and received packets, drawn from some ambient 
message space fJ, can be gathered into the rows of a transmit- 
ted matrix X and a received matrix Y . Error packets injected 
into the network can be described by the rows of a noise 
matrix Z. Due to the nature of linear network coding, the 
linear transformation of transmitted packets X and the linear 
propagation of error packets Z can be modelled as an additive- 
multiplicative matrix channel (AMMC), defined via 



Y = AX + BZ 



(1) 



for some transfer matrices A, B. One typically assumes that 
A, B, and Z are random matrices (drawn according to some 
distributions) and independent of X. This type of stochastic 
model is appropriate in situations where random network 
coding is performed and the noise Z arises due to decoding 
errors, rather than from the malicious actions of an adversary. 
When the ambient space fi is a vector space over a 
finite field, tight capacity bounds and simple, asymptotically 
capacity-achieving, coding schemes are developed in jT], 
under certain distributions of A, B, and Z. Similar work along 
this line can be found, e.g., in |l2l-||5]. Prior work on matrix 
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channels for linear network coding has mainly focused on the 
finite-field case. 

In this paper, we consider the case when the ambient space 
is an arbitrary finite module over some finite chain ring, 
generalizing the work of [l] from finite fields to finite chain 
rings. The motivation for considering this generalization arises 
from nested-lattice physical-layer network coding ll6l- lfT0l . in 
which the ambient space 51 is a lattice quotient A/ A', where 
A is some complex integer lattice and A' is a sublattice. 
Algebraically, such an ambient space is a finite module over 
some principal ideal domain (PID) T of the form 

n^T/{di) xT/{d2) x...xT/(d„), 

where di,d2, ■ ■ ■ , dm G T are nonzero non-unit elements 
satisfying the divisibility conditions di | ^2 | • • • | dm- The 
PID T is often chosen as a discrete subring of C, such as 
the integers Z, the Gaussian integers Z[i], or the Eisenstein 
integers Z[e^'^'^^^]. In many useful lattice constructions (see, 
e.g., mi ), the ambient space il takes the special form 

n = T/{p'')x---xT/{p'-), (2) 

where p is a prime in T, and ti, . . . ,tm are positive integers 
satisfying ii < • • • < i^, ||8]. In this case, Q, is isomorphic to 
a finite module over the finite chain ring r/(p*™). 

It is thus of interest to consider matrix channels over finite 
chain rings. As in |r|, we gather insight by first studying two 
variations: the noise-free multiplicative matrix channel (MMC) 
Y = AX, and the multiplication-free additive matrix channel 
(AMC) Y = X + BZ. 

The essential step in handling the MMC over finite fields 
is to use the so-called principal Schubert cell that consists 
of subspaces that occur as the row space of all matrices 
having a certain reduced row echelon form (RREF) ([T]. The 
extension of these concepts is not straightforward; the main 
problem is the presence of zero divisors in the ring. For 
example, one existing extension of the RREF is called the row 
canonical form [12], which, however, only applies to square 
matrices over Galois rings (a special case of finite chain rings). 
Another extension is called the Howell form lfT3l . which, 
however, fails to give a concise representation (in terms of 
the number of nonzero rows). To overcome these difficulties, 
we introduce a generalized notion of row canonical form for 
matrices of any size with elements from a finite chain ring 
and which contains the minimum number of nonzero rows. 
A generalization of the Gauss-Jordan elimination procedure 
allows one to convert any given matrix, using elementary row 
operations, to a matrix in row canonical form having the same 
row module. Based on this notion, we then give a natural 



generalization of the principal Schubert cell, namely, the set 
of modules that occur as the row module of all matrices having 
a certain row canonical form; furthermore, we show how this 
set of modules can be constructed efficiently. Equipped with 
these generalizations, we are able to obtain results parallel 
those of [H. 

The key step in handling the AMC over finite fields is 
to use a so-called error-trapping strategy [1]. In finite-field 
matrix channels, the size of the error-traps is determined by 
the number of matrices of a given rank t. The rank t may 
be regarded as a measure of "noise level" of the matrix BZ. 
For matrices over finite chain rings, the concept of "rank" is 
more subtle, and must be suitably generalized. We first show 
that how the concept of "shape" — the appropriate chain-ring- 
theoretic generalization of dimension — can be used to indicate 
the noise level. We then derive an enumeration result on the 
number of matrices of a given shape, which enables us to 
determine the size of the error-traps needed in our coding 
scheme. This enumeration result generalizes that of llT4l from 
square matrices to general matrices and from Galois rings to 
finite chain rings, and plays a crucial role in obtaining the 
capacity and the capacity bounds in the finite-chain-ring case. 

Building upon the generalizations for the two special cases, 
we derive tight capacity bounds and simple, polynomial- 
complexity, asymptotically capacity-achieving coding schemes 
for AMMC models related to ^. 

The remainder of this paper is organized as follows. Sec- 
tion [n] reviews some basic facts about rings, modules and 
matrices over finite chain rings. Section |III] introduces the 
row canonical form. Section |IV] presents several enumeration 
results and construction methods for matrices under row 
constraints. The new results developed in Sections Hill and HVl 
serve as building blocks to the subsequent sections. Section IV] 
formalizes the basic matrix channels that are studied in this 
paper. The three variations (MMC, AMC, and AMMC) are 
addressed in Sections [Vl] IVIII and IVIIII respectively, where 
capacity and coding results are presented. Finally, Section |IX] 
presents our conclusions. 

II. Preliminaries 

In this section, we present some basic results for finite chain 
rings, and modules and matrices over finite chain rings. This 
section establishes notation and the results that will be used 
later; nevertheless, this material is standard; see e.g., lfT2ll . 
lUSj-IIll for more details. 

A. Rings and Ideals 

We begin with some common definitions and notations for 
rings. All rings in this paper will be commutative with identity 
1 7^ 0. Let i? be a ring. We will let R* denote the nonzero 
elements of R, i.e., R* = R\{0}. An element a in i? is called 
a unit if a6 = 1 for some b E R. We will let U{R) denote the 
units in R. Two elements a,b E R aie said to be associates if 
a = ub for some u E U{R). Associatedness is an equivalence 
relation on R. 

Suppose a,b E R. The element a divides b, written a \ b, 
if ac — b for some c E R. Let d E R* he a. nonzero element 



in R. Two elements a, b are said to be congruent modulo d 
if d divides a — b. Congruence modulo d is an equivalence 
relation on R. A set containing exactly one element from 
each equivalence class is called a complete set of residues 
with respect to d, and is denoted by TZ{R, d). Note that the 
difference a — b between distinct elements a,b E TZ{R, d), 
a y^ b, can never be a multiple of d. 

An element a of R* is a called a zero-divisor if a6 = 
for some b E R*. If R contains no zero-divisors, then R is 
an integral domain. If R is finite and an integral domain, then 
R is, in fact, a finite field. This latter case is not of central 
interest in this paper; almost all of the rings considered here 
will have zero divisors. 

Example 1: Let R = Zg = {0, . . . , 7}, under integer addi- 
tion and multiplication modulo 8. Then U{Zg,) = {1, 3, 5, 7}. 
There are four equivalent classes induced by congruence mod- 
ulo 4, namely, {0,4}, {1,5}, {2,6}, and {3,7}. An example 
of a complete set of residues with respect to the element 4 
in TZ{Zs) is 7?.(Z8, 4) = {0, 1, 2, 3}. The zero-divisors of Zg 
form the set {2, 4, 6}. ■ 

A nonempty subset I of R that is closed under subtraction, 
i.e., a.b E I implies a — b E I, and closed under inside-outside 
multiplication, i.e., a E I and r E R implies ar E I, is called 
an ideal of R.lf A = {ai, . . . , a^} is a finite nonempty subset 
of R, we will use (ai, . . . , am) to denote the ideal generated 
by A, i.e.. 



(ai, 



i) = {aici H V amCm'- Ci, . . .,€,„ E R}. 



An ideal / of i? is said to be principal if / is generated by a 
single element in /, i.e., / = (a) for some a G /. A ring R 
is called a principal ideal ring (PIR) if every ideal / of i? is 
principal. If i? is a PIR and also an integral domain, then R 
is called a principal ideal domain (PID). 

An ideal A^ is said to be maximal if N ^ R and the 
only ideals containing N ait N and R (in other words, N 
is "maximal" with respect to set inclusion among all proper 
ideals). If iV is a maximal ideal, then the quotient R/N is 
a field, called a residue field. A ring with a unique maximal 
ideal is said to be local. 

Example 2: The ideals of Zg are {0} = (0), {0,4} = 
(4), {0,2,4,6} = (2), and R = (1). Thus, Zg is a PIR, and 
has a unique maximal ideal (2). The residue field Zg/(2} is 
isomorphic to the finite field F2 of two elements. ■ 

B. Finite Chain Rings 

A ring R is called a chain ring if the ideals of R satisfy 
a containment condition: for any two ideals /, J of R, either 
/ C J or J C /. If i? is a chain ring with finitely many 
elements, then R is called a finite chain ring. Clearly, a finite 
chain ring has a unique maximal ideal, and hence is local. It 
is known lfT2l that a finite ring is a chain ring if and only 
if it is a local PIR; thus, in a finite chain ring, all ideals are 
principal. Examples of finite chain rings include Zp>. (the ring 
of integers modulo p" where p is a prime) and Galois rings. 

Let i? be a finite chain ring, and let tt € i? be any generator 
of the maximal ideal of R. Then R/{tt) is the residue field 



of R. It can be shown (see, e.g., lfT2l ') that every ideal I of 
R, including the zero ideal (0), is generated by a power of vr, 
i.e., / = (tt') for some / > 0. It follows that tt is nilpotent; we 
denote by s the nilpotency index of tt, i.e., the smallest positive 
integer such that tt^ — 0. There are, then, exactly s + 1 distinct 
ideals of R, namely, R = (tt"), (tt^), . . . , (tt'*) = {0} which 
form a chain (with respect to set inclusion): 

R = (tt") D (ttI) D • • • D (vr^-i) D (vr^) = {0}. 

Thus, s is often called the chain length of i?. We refer to R 
as a (g, s) chain ring, if R has a residue field of size q and a 
chain length of s. 

Example 3: The ideals of Zg form a chain with respect to 
set inclusion: 

i?=(l)D(2)D(4)D(0) = {0}. 

Thus, Zg is a finite chain ring with chain length s = 3. Since 
the residue field Z8/(2) is isomorphic to F2, Zg is a (2,3) 
chain ring. ■ 

Now let Tl{R, tt) C _R be a complete set of residues 
with respect to tt and, without loss of generality, assume 
that G TZ{R,tt). Every element a E R then has a unique 
representation, called the ir-adic decomposition of a (with 
respect to Tl{R, tt)), in the form 



ao + aiTT + • • • + tts-iTT 



where oq, . . . , as-i G TZ{R, tt). It follows from the uniqueness 
of (O that the size of R is q^, i.e., the number of elements 
in a (g, s) chain ring is q^ . Thus, like a finite field, a finite 
chain ring has a cardinality that is an integer power of a prime 
number 

The degree of a nonzero element aq + oitt 
as-iTT*"^ G i?*, denoted by deg(a), is defined as the least 
index j for which aj 7^ 0. By convention, the degree of is 
defined as s. All elements of the same degree are associates 
in R. Further, a divides h if and only if deg(a) < deg(6). 
Finally, deg(a + &) > min{deg(a),deg(&)}, i.e., adding two 
elements cannot result in an element of lower degree. 

Example 4: Let 7e(Zg,2) = {0,1}. The 2-adic decompo- 
sition of 5 G Zg is 5 = 1 + ■ 2 + 1 ■ 2^. The elements 
in Zg of degree (respectively, 1, 2, and 3) are {1,3,5,7} 
(respectively, {2, 6}, {4}, and {0}). ■ 

Finally, we present two methods for constructing finite chain 
rings. 

If R is itself a (g, s) chain ring with maximal ideal (vr), 
then the quotient R/{tt^) (0 < / < s) is a {q,l) chain ring. 
This method constructs new finite chain rings from existing 
ones. 

If T is a PID, and p is a prime in T, then T/{p) is a field, 
since (p) is a maximal ideal of T. Let q be the size of T/{p) 
and suppose that q is finite. Then the quotient T/ {p^) is a (q, /) 
{I > 0) chain ring. This method constructs finite chain rings 
from PIDs. 



C. Modules over Finite Chain Rings 

A module is to a ring as a vector space is to a field. More 
formally, an _R-module A/ is an abelian group (A/, +) together 
with an action of R on M satisfying the following conditions 
for all TO, n G M and for all a, & G -R: 

1) Ito, = m and {ab)m — a{bm) 

2) (a + b)m — am + bm 

3) a{m + n) — am + an 
When R is a finite chain ring, an i?-module is always 

isomorphic to a direct product of various ideals of R; this 
structure can be described by a "shape." An s-shape /i = 
(/Ui, /i2, ■ . ■ , lis) is simply a sequence of non-decreasing non- 
negative integers, i.e., < /ii < /i2 < • ■ • < Ms- We denote by 
\fj,\ the sum of its components, i.e., \fi\ — ^21=1 l^i- ^'^^ ^^^^^ 
notational convenience, we define the "zeroth component" of 
a shape as /ip — 0. 

An s-shape k = {ki, . . . ,Ks) is said to be a subshape 
of /i = (/ii, . . . ,/is), written k ^ /i, if k^ < jii for all 
i = l,...,s. Thus, for example, (1,1,3) ^ (2,4,4). The 
number of subshapes of the s-shape (m, . . . , m) is given 
by (™j^'*), which implies that the number of subshapes of 
M = (^1, . . . , /Xs) is upper-bounded by (^%"^*). 

Two s-shapes can be added together to form a new s- 
shape simply by adding componentwise. Thus, for example, 
(1,1,3) + (2,4,4) = (3,5,7). Also, for a shape ^jl = 
(/ii, . . . ,/is) and a positive integer to. we define ji/m, = 
(/ii/m, . . . ,/is/TO,) (which is an s-tuple, but not necessarily 
(■') a shape). For convenience, we will sometimes identify the 
integer t with the s-shape (t, . . . , i). Thus, for example, /i ^ t 
means Hi < t for all i, k = t means kj = t for all i, and 
fi — t — (/ii — t, . . . , fis — t), assuming t < ji. 

Let i? be a (g, s) chain ring with maximal ideal (tt). For 
any s-shape /i, we define the _R-module _R'' as 

+ ai-K + ■■■ + i?^' A (|;l^ X • . • X (1) X (tt) X . • • X (tt) X • • • 

/J2— Ml 

-S-1 



(tt'^-I) X •••X (tt'^-I). (4) 



Since a positive integer t is identified with the shape (t, . . . , t), 
it is indeed true that i?* denotes the t-fold Cartesian product 
of R with itself. 

The module _R^ can be viewed as a collection of ^s- 
tuples whose components are drawn from R subject to certain 
constraints imposed by /i. Specifically, while the first /ii 
components can be any element of R, the next /i2 — /ii 
components must be multiples of tt, and so on. Since each 
ideal (tt*) in (HI contains q^^^ elements (0 < i < s), it follows 
that the size of i?^ is |i?^| = gl^L 

Example 5: Let R = Zg, and let ^ = (2,4,4). Then 

i?'' = (1) X (1) X (2) X (2) . 



Note that the first two components of i?'^ can each be chosen 
in 2^ ways, while the last two components can each be chosen 
in only 2^ ways. Hence, the size of i?^ is 2^". ■ 



For every s-shape ^, R^ is a finite i?-module. Conversely, 
the following theorem establishes that every finite i?-module 
is isomorphic to R^ for some unique s-shape /i. 

Theorem 1: |18, Theorem 2.2] For any finite i?-module 
M over a {q, s) chain ring R, there is a unique s-shape /i such 
that M ^ R^'. 

We call the unique shape /i given in Theorem [T]f/ie shape of 
M, and write /i = shape MJ^ It is known |18| that if M' is a 
submodule of M, then shape A/' ^ shape M, i.e., the shape 
of a submodule is a subshape of the module. It is also known 
ifTSl that the number of submodules of R^ whose shape is n 
is given by 



H' 



where 



= 9 



,{P-i-Hi)K,i 



(5) 



n 

i=0 



is the Gaussian coefficient. In particular, when the chain length 
s = 1, i? becomes the finite field Fg of q elements, and 
[^] becomes [^^] , which is the number of ki -dimensional 
subspaces of F''^ . 

D. Matrices over Finite Chain Rings 

We turn now to matrices over finite chain rings. Let i? be a 
{q, s) chain ring with maximal ideal (tt). The set of all n x m 
matrices with entries from R will be denoted by i?"^™. If 
A £ i?"^™, we denote by A[i,j] the entry of A in the zth 
row and jth column, where 1 < i < n and I < j < rn. We 
will let A[ii:i2, ji:J2] denote the submatrix of A formed by 
rows zi to «2 and by columns ji to J2, where 1 < ii < 12 < n 
and 1 < ji < J2 < rn. Finally, we will let A[i,:] denote the 
ith row of A and ^[:, j] denote the jth column A. 

A square matrix [/ e i?"^" is invertible if UV — VU = /„ 
for some V G i?"^", where /„ denotes the n x n identity 
matrix. The set of invertible matrices in i?"^", denoted as 
GL„(i?), forms a group — the so-called general linear group — 
under matrix multiplication. 

Two matrices A,B^ j^nxm ^.g ^^^^ j.^ ^g left-equivalent 
if there exists a matrix U G GL„(_R) such that UA — B. 
Two matrices A, _B G j^nxm ^j.g ^^j^ ^^ |^g gquiyalent if there 
exist matrices U G GL„(i?) and V G GL„i(i?) such that 
UAV = B. 

A matrix D G i?"^™ is called a diagonal matrix if I?[i, j] — 
whenever i ^ j. A diagonal matrix D, which need not 
be square, can be written as 13 = diag((ii, . . . , d^), where 
r = min{n, to}, and di — D[i, i] for ? = 1, . . . , r. 

Let A G R^^"\ A diagonal matrix D = diag((ii, . . . ,dr) G 
j^nxm ^^ _ niin{n,TO,}) is called a Smith normal form of A, 
if D is equivalent to A and di | (i2 | ■ • • | dr in i?. It is known 
IIT6I that every matrix over a PIR (in particular, a finite chain 
ring) has a Smith normal form whose diagonal entries are 

'Some authors (like Honold et al. 1 181 ) use a different convention and 
define the shape of an i?-module to be the conjugate (in the integer-partition- 
theoretic sense) of the shape as defined in this paper. 



unique up to equivalence of associates. In this paper, we shall 
require the diagonal entries di, . . . ,dr in the Smith normal 
form D to be powers of tt, i.e., 

(di,...,4) = (V\...,V'-), 

where < li < . . . < Ir < s since di \ d2 \ ■ ■ ■ \ dr- With 
this constraint, once vr is fixed, every matrix A G i?"^'" has 
a unique Smith normal form. 

Example 6: Consider the two matrices 



A = 



4 6 


2 1' 
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6 1 
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over Zg. It is easy to check that 
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= usv. 



Since U and V are invertible, S is equivalent to A. Since the 
diagonal entries of S satisfy 1 | 2 | 4 | in Zg, S is the Smith 
normal form of A. ■ 

For any A G i?"^™, we denote by rowyl and col^ the 
row span and column span of A, respectively. By using the 
Smith normal form, it is easy to see that the row span row A 
is isomorphic, as an i?-module, to the column span colyl. 

Note that two matrices A,Be j^nxm ^j.g left-equivalent 
if and only if row A = towB, i.e., left-equivalent matrices 
have identical row spans. On the other hand, two matrices 
A, i? G R"-'""'- are equivalent if and only if row A = row B, 
i.e., equivalent matrices have isomorphic row spans. 

The shape of a matrix A is defined as the shape of the row 
span of A, i.e., 

shape A = shape (row A). 

Clearly, shape A — shape(colyl). Moreover, shape A = /i if 
and only if the Smith normal form of A is given by 



diag(l,...,l,7r, ...,7r, 



.,7r 



.,^^-1,0,, 



,0), 



where r = min{n,TO.}. In particular, a matrix U G i?"^" is 
invertible if and only if shape U — {n, . . . ,ti). 

Example 7: Since D = diag(l, 2, 4, 0) is the Smith normal 
form of A in Example |6] shape A ~ (1,2, 3). ■ 

As one might expect, matrix shape has a number of prop- 
erties similar to matrix rank. 

Proposition 1: Let A G R"""'"- and B e R'"'"'. Then 

1) shaped = shape A^, where A^ is the transpose of A. 

2) For any P G GL„(i?),Q G GL,„(i?), shaped = 
shape PAQ. 

3) shape AS ^ shaped, shape AS ^ shape B. 

4) For any submatrix C of A, shape C :< shape A. 



Proof: 1) Since row A = col A, we have row A = 
rowA"^. Hence, shape A = shape A-^. 2) Since A is equiv- 
alent to PAQ for any invertible P and Q, shape A — 
shape PAQ. 3) Since row AS is a submodule of rowB, we 
have shape A_B ^ shape B. Similarly, since col AB is a sub- 
module of col A, we have shape AB :< shape A. 4) Finally, 
note that any submatrix C of A is equal to E1AE2 for some 
El e R^^^ (selecting k rows) and E2 £ R"^^'- (selecting I 
columns). Hence, shape C — sYv&Yie E1AE2 di shape A. ■ 

The shape of A is closely related to the McCoy rank of A 
lUSl, ri6|. Specifically, the (McCoy) rank of A is equal to the 
first component of the shape of A, i.e., rank A — /ii, where 
/^ = shape A. 

A mati-ix A e j^nxm j^ called full rank if rankvl — 
min{n,m}. A matrix A £ Jinxm j^ called full row rank if 
rankA = n (which requires n < m). The number of full- 
row-rank matrices in R"^^™ is qS"™ nr=o (1 — 9*~™)- Clearly, 
A e /j"x™ is full row rank if and only if shape A = n. A 
matrix is full column rank if its transpose is full row rank. 

A famous result, due to N. McCoy, connects the rank of a 
matrix A with the system of equations AX = 0. 

Theorem 2: lUg] Theorem 5.3] Let A e i?"^™. The 
system of equations AX = has a nontrivial solution if and 
only if rank A < ra. 

III. Row Canonical Form 

The row canonical form presented here slightly generalizes 
the row canonical form in |12, Exercise XVI. 7] from Galois 
rings to arbitrary finite chain rings and from square matrices to 
arbitrary matrices. As we will see, this new row canonical form 
is a natural analogue, for matrices over finite chain rings, of 
the familiar reduced row echelon form for matrices over fields. 
In addition to describing an algorithm for reducing an arbitrary 
matrix to row canonical form, we provide an elementary proof 
for the uniqueness of this form. 

Throughout this section, i? is a (g, s) chain ring with 
maximal ideal (tt). We fix a complete set of residues TZ{R^ vr) 
(including 0), i.e., a representation of the residue field R/ {tt), 
and, for 1 < / < s, we choose the complete set of residues 
for tt' as 

n{R, tt') - < ^ a.y : ao, . . . , a^-i G 7^(i?, tt) i . 
Finally, we set TZ{R,tt^) = {0}. 

A. Definitions 

We start with a few definitions. 

Let A be matrix with entries from R. The ith row of A 
is said to occur above the (i')th row of A (or the (i')th row 
occurs below the ith row) if i < i'. Similarly the jth column 
of A is said to occur earlier than the (j')th column (or the 
(j')th column occurs later than the jth column) if j < j' . 
This terminology extends to the entries of A: A[i,j] is above 
A[i' ,j'] if i < i' and A[i,j] is earlier than A[i',j'] if j < j'. 
If P is some property obeyed by at least one of the entries in 
the ith row of A, then the first entry in row i with property P 
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occurs earlier than every other entry in row i having property 
P. 

The pivot of a nonzero row of a matrix is the first entry 
among the entries having least degree in that row. For example, 
6 and 2 are the entries of least degree in the row [0 4 6 2] 
over Zg, and 6 occurs earlier Thus, 6 is the pivot of the row 
[0462]. Note that the pivot of a row is not necessarily the 
first nonzero entry of the row. 

Definition 1: A matrix A is in row canonical form if it 
satisfies the following conditions. 

1) Nonzero rows of A are above any zero rows. 

2) If A has two pivots of the same degree, the one that 
occurs earlier is above the one that occurs later. If A 
has two pivots of different degree, the one with smaller 
degree is above the one with larger degree. 

3) Every pivot is of the form tt' for some I G {0, . . . , s — 1}. 

4) For every pivot (say tt'), all entries below and in the 
same column as the pivot are zero, and all entries above 
and in the same column as the pivot are elements of 
n{R,Tr^). 

Example 8: Consider the matrix 



A 



over Zg with tt = 2 and Z8/(2) = {0,1}, in which the pivots 
have been identified with an overline. Clearly, A satisfies all 
of the conditions to be in row canonical form. ■ 

The following facts follow immediately from the definition 
of row canonical form. 

Proposition 2: Let A £ r^x™- be a matrix in row canonical 
form, let p^ be the pivot of the fcth row, let c^ be the index of 
the column containing p^. (If the fcth row is zero, let p^ — 

and Cfc = 0.) Let dk = deg{pk), and let w = (wi, . . . , Wm) 
be an arbitrary element of row A. 

1) Any column of A contains at most one pivot. 

2) If A has more than one row, deleting a row of A results 
in a matrix also in row canonical form. 

3) i > k implies deg{A[i,j]) > dk- 

4) (i > k and j < Ck) or (i > k and j < Ck) implies 
deg{A[i,j]) >dk. 

5) pi divides wi,W2,- ■■ ,ujjn. 

6) j < ci implies deg(wj) > di. 

Proof: 1) The presence of a pivot p in a column rules 
out the possibility of another pivot in the same column and 
below p, since all entries in the same column below p must 
be zero and hence cannot be pivots. 2) Deleting a row of A 
does not influence the value or the position of the pivots in 
the other rows; thus it easy to verify that the modified matrix 
satisfies the four conditions required for a matrix to be in row 
canonical form. 3) By definition pk has degree smaller than 
or equal to that of any element in its row. If A contained 
an element in a row below row k of degree smaller than dk, 



then the pivot of that row would have degree smaller than dk, 
contradicting the property that pivots of smaller degree must 
occur above pivots of larger degree. 4) By definition pk is the 
earliest element having minimum degree in row k, so every 
element in row k occurring earlier than pk has degree strictly 
larger than dk- We know from 3) that A contains no element 
in a row below k of degree smaller than dk- If such a row 
contains an element of degree equal to dk, then the pivot of 
that row must occur later than pk, which implies that every 
element occurring in that row occurring in column Ck or earlier 
has degree strictly larger than dk- 5) Consider Wj- From 3) 
we know that pi divides every element of A; in particular, pi 
divides every element of column j of A- Since Wj is a linear 
combination of these elements, it must be that pi divides Wj . 
6) If j < c\, we know from 4) that every element in column 
j of A has degree strictly greater than d\ and so does every 
linear combination of these elements, in particular Wj- ■ 

For any A E i?"^™, we say a matrix B G Ji^xm j^ ^ ^^^ 
canonical form of A, if (i) B is in row canonical form, and (ii) 
B is left-equivalent to A- We will show that any A G _R"^™ 
has a unique row canonical form. For this reason, we denote 
by RCF(j4) the row canonical form of A- 

B. Existence 

We will demonstrate the existence of a row canonical 
form for any matrix A by presenting a simple algorithm that 
performs elementary row operations to reduce A into row 
canonical form. Here, the allowable elementary row operations 
(over K) are: 

• Interchange two rows. 

• Add a multiple of one row to another 

• Multiply a row by a unit in R. 

Each of these operations is invertible, and so a matrix obtained 
from A by any sequence of these operations will have the same 
row span as A- 

The algorithm proceeds in a series of steps. In the fcth step, 
the algorithm selects the fcth pivot, moves it to the fcth row, and 
uses elementary row operations to reduce into row canonical 
form the submatrix consisting of the top k rows. The pivot 
selection procedure operates on any given set of rows. If the 
rows are all zero, the procedure should return with the result 
that no pivot can be found. Otherwise, among all entries of 
least degree in the given rows, an entry must be chosen that 
occurs as early as possible. This entry must certainly be the 
pivot of its row. The procedure should return the row and 
column index of the selected element. 

Now we are ready to describe the algorithm in detail. In 
step A; = 1, apply pivot selection to all of the rows of A. If 
no pivot can be found, then A is a zero matrix, and is already 
in row canonical form. Otherwise, we call this pivot the first 
pivot and place it in the first row by an interchange of rows (if 
necessary). If this pivot is not of the form tt' (/ = 0, . . . , s — 1), 
we multiply the first row by a suitable unit so that the first 
pivot is a power of tt. Note that nonzero entries in the same 
column below the first pivot have degrees no less than the 
pivot, which means that they are all multiples of the first pivot. 
By a sequence of elementary row operations, these entries can 



be cancelled, so that we arrive at a matrix, say Ai, in which the 
first row is in row canonical form and all entries in the same 
column below the first pivot are zero. We can now increment 
k and proceed to the next step. 

For fc > 2, we apply pivot selection to the rows of Ak-i, 
excluding the first fc — 1 rows. If no pivot can be found, then 
the remaining rows are all zero and Ak-i is in row canonical 
form. Otherwise we call this pivot the fcth pivot and place it 
in the fcth row by an exchange of rows (if necessary). As in 
the first step, if this pivot is not an integer power of tt, we 
multiply the fcth row by a suitable unit so that the fcth pivot 
is a power of tt, say tt'. Nonzero entries in the same column 
below the fcth pivot can be cancelled using elementary row 
operations. A nonzero entry, say a, in the same column above 
the fcth pivot has 7r-adic decomposition 



a = flo + • • • + fls-lTT 



ao 



aj_i7r + TT (aj + • • • + as_i7r 



s-;-i^ 



Thus by subtracting (a; + • • • + aj-iTr"^ ) times the fcth 
row from the row containing a, we change a to gq + ■ ■ ■ + 
a;_i7r'~^ G TZ{R, tt'), without affecting the pivot of that row. 
Reducing all nonzero entries in the same column as the fcth 
pivot in this way, we arrive at a matrix, say Ak, in which the 
top fc rows are in row canonical form and all entries in the 
same column below the first, second, . . . , fcth pivots are zero. 

The above algorithm stops when no more pivots can be 
found. Note that, at the end of the fcth step, the matrix Ak is 
left-equivalent to A and the submatrix formed by the top fc 
rows of Ak is in row canonical form. It follows that the final 
matrix must be in row canonical form. 

Therefore, we have the following result. 

Proposition 3: For any A G i?"^™, the algorithm de- 
scribed above computes a row canonical form of A- 

A simple count shows that this algorithm requires 

0{nm niin{n, m}) 

basic operations over R- 

Example 9: Consider the matrix 



A 



over Zg. There are three Is in the last column of A, namely, 
A [1,4], A [3, 4] and A [4, 4], which are the elements of least 
degree in A- We can choose any of them as the first pivot. 
Here, we choose A[l, 4] (indicated by an overline). After some 
elementary row operations, we can make the entries below the 
pivot zero to obtain 
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Now consider the submatrix formed by omitting the first 
row of Ai- There are four entries of least degree, namely. 
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Ai[3,l] = 6, Ai[3,2] = 6, Ai[4,l] = 6, and Ai[4,2] = 2, 
among which ^i[3,l] and yli[4, 1] are valid choices for 
the second pivot. Here, we choose Ai[3, 1] (indicated by an 
overUne). We interchange the second row and third row of Ai, 
and then muhiply the new second row by 3, obtaining 



A[ 



By some elementary row operations, we can make the entries 
below the second pivot zero. After that, we subtract 2 times 
the second row from the first row, obtaining 



A,= 



Clearly, the submatrix formed by the top two rows of A2 is in 
row canonical form. Next, consider the submatrix formed by 
omitting the top two rows of A2. We choose the entry A2[3, 2] 
(indicated by an overline) as the third pivot. We subtract the 
third row from the fourth row and obtain 



A. 



Clearly, the submatrix formed by the top three rows of A3 is 
in row canonical form (with all the pivots indicated). Since no 
more pivots can be found, our algorithm outputs A3, which is 
indeed in row canonical form. ■ 



C. Uniqueness 

We will prove the uniqueness of the row canonical form 
here. 
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Proposition 4: For any A e R"^ 
of A is unique. 



\ the row canonical form 



Proof: If A is the zero matrix, then its row canonical 
form must also be the zero matrix, which is therefore unique. 
Thus let us assume that A is nonzero. 

We will proceed by induction on n. For n — 1, the proof 
is obvious. Thus suppose that n > 1, and let B and C be 
two row canonical forms of A. Clearly, row B = row C, and 
each row of B and C are elements of row A. Let B[l,ji] and 
C[l, J2] be the pivots in the first row of B and C, respectively. 
From Proposition |2]-5 we have that B[l,ji] \ C[1,J2] and 
C[1,J2] I i?[l,ji]; thus i?[l,ji] and C[1,J2] are associates. 
However, since pivot elements must take the form tt' for some 
I, we conclude that B[l,ji] = C[1,J2]- Suppose ji < J2- 
By Proposition |2]-6 we have deg(_B[l, ji]) > deg{C[l,J2]), 
contradicting the fact that B[l,ji] = C[1,J2]- A similar 
contradiction arises if ji > j2- We conclude that ji — J2, 
i.e., both B and C must have exactly the same pivot element 
in exactly the same position in their first row. 



Now let ji = J2 = j- Consider the submodule of rowyl 
in which every element has zero in its jth component. Every 
element a of this submodule is a linear combination 

n 

a = y^^biB[i,:]; 

for some choice of coefficients bi,...,bn- However, since 
flj = 0, and B[i,j] = for i > 2, we must have biB[l,j] = 0. 
Since B[l,j] is the pivot element of the first row of B, it 
divides every element of that row; thus if biB[l,j] = 0, 
then biB[l,:] = 0, i.e., the first row can only contribute 
to a. This means that the given submodule is equal to 
row_B[2:n, l:m]. Similarly, the given submodule is also equal 
to rowC[2:n, l:m]. By Proposition |2]-2, both B[2:n, 1:to] and 
C[2:n, 1:to] are in row canonical form. Thus by induction, we 
have B[2:n,l:'m] — C[2:n,l:m]. This implies that B and C 
can differ in their first row only. 

Let us assume that _B[1, :] j^ C[l, :], i.e., that the first rows 
of B and C are not equal, so that A — {61, . . . ,dm) = 
B[l, :] — C[l, :] is nonzero. Since A is an element of row A 
with zero in its jth component, we have A G row i?[2:n, l:m], 
from which it follows that 



A = ^c,B[z,:] 



1=2 

for some C2, . . . , c„ G R. If B[2:n, l:m] is the zero matrix, 
then A = 0, which is a contradiction. Otherwise, let B[2,J3] 
be the pivot of i?[2, :]. Note, on the one hand, that _B[i, J3] = 
for all i > 2; thus 6j^ = C2i?[2, J3], i.e., 5j^ must be a multiple 
of B[2,J3]. On the other hand, because _B[2, J3] and C[2,j3] 
are (identical) pivots, B[l,j3], C'[1,J3] e TZ{R, B[2,J3]). If 
B[1,J3] and C[l,j3] are distinct, their difference, 5j^, cannot 
be a multiple of i?[2,j3]. We conclude that 6j^ = 0, i.e., 
B[1,J3\ and C[1,J3] are not distinct. Since _B[2,J3] is the 
pivot of i?[2,:] it divides every element of -B[2, :]; thus if 
C2i?[2, J3] — 0, then C2i?[2, :] = 0. Continuing this argument, 
we have CiB[i, :] = for all i > 2. Therefore, we have A = 0, 
which is a contradiction. This establishes uniqueness. ■ 



D. Applications 

We will provide two immediate applications of row canon- 
ical forms. 

1) Computing the Shape of a Matrix: Previously, the shape 
of a matrix was obtained from the Smith normal form. Here, 
we will connect the row canonical form with the Smith normal 
form so that the shape can be computed from the row canonical 
form as well. 

The connection is as follows. Let B be the row canonical 
form of A G _R"^™ with k nonzero rows. Let pi be the pivot in 
the ith row of B, where i G {!,..., k}. Let r = min{n, m}. 
Clearly, k < r. Then the Smith normal form of A is given by 

diag(pi,...,Pfe,0,...,0)Gi?"^'", 

r—k 

from which the shape of A is readily available. 



This connection has another interpretation. For any a ^ R, 
we define the modulo-7r* (i — 1 , . . . , s) operation as 



a mod TT* 



(tt'). 



This operation extends to matrices in an element-wise fashion 
so that B mod tt* gives a matrix over i?/(7r'). Let fi — 
shape A. Then /ij is the number of nonzero rows in B mod 
TT* (over _R/(7r*)), for i = 1, . . . , s. 

Example 10: Consider the following two matrices 
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over Zg. Since B is the row canonical form of A, and the 
number of nonzero rows in B mod 2 (respectively, mod 4, 
mod 8) is 1 (respectively, 2, 3), we conclude that shape A — 
(1,2,3). ■ 

As a simple consequence, the rank of A is the number of 
nonzero rows in B mod tt. This implies that A is full rank 
(over R) if and only if A mod vr is full rank (over R/{tt)). 

2) Rank Decompositions: We introduce a decomposition 
related to the row canonical form, which can be viewed as a 
natural extension of the well-known rank decomposition. 

Proposition 5: Let B be the row canonical form of A S 
i?"^™. Let B be the submatrix of B consisting of only 
nonzero rows. Then A can be decomposed as a product 
PiB of some full-column-rank matrix Pi and the matrix B. 
Moreover, the number of Pi producing such a decomposition 



IS q 



" Ei=i iif^i-t 



where k — shape A. 



Proof: Since B is the row canonical form of A, A = 
PB for some invertible matrix P £ GL„(i?). Since k = 
shape A — shape i?, B has Kg nonzero rows, and B G 
^K,xm Let P = [Pi Pa], where Pi G P"^'*^^ and 
P2 G P"x(«-«=). Then we have 



A = PB= [Pi P2] 



PiP. 



Since P is invertible, shape P = {n,. . . ,n). In particular, P 
mod TT is full rank over R/ {tt). It follows that Pi mod vr is full 
column rank over R/{tt). Hence, shape Pi = (ks,...,Ks), 
which implies that Pi is full column rank over R. 

Next, we count the number of such decompositions. Con- 
sider the matrix equation XB = PiB, in unknown X. Clearly, 
the number of decompositions of A is equal to the number 
of solutions to this matrix equation. Let B[i,ji] be the pivot 
of the ith row of B, for all i = 1,...,Ks. Then B[i,ji] 
divides the ith row of B. It follows that B = DB', where 
D = diag I P[l, ji], . . . ,B[ks,JkJ ), and the ith row of B' 

is equal to the ith row of B divided by B[i,ji]. Clearly, 
B'[i,ji] = 1 for all i = l,...,^^. Since ji,...,Jk, are 
all distinct, B' mod vr is full row rank over R/ {tt). Hence, 
shape P' = (ks,...,Ks), which implies that B' is full row 
rank over R. 

Since B' is full row rank, by TheoremH {XD-PiD)B' = 
if and only if XD- PiD ^ 0. Hence, XB ^ PiB if and 



Xu- 




X, 



< Ml > < /'2 ► 

Fig. 1 . Illustration of a vr-adic decomposition for 




3 and fi = (4, 6, 



only if XD = PiD. Thus, it suffices to count the number of 
solutions to XD == PiD. Note that XD = PiD is equivalent 
to the following system of equations 



X[i,k]B[k,jk] ^ Pi[i,k]B[k,jk],i = 1, 



, n, fc = 1, , 



(6) 
Suppose that B[k,jk] — tt'-'' for some < fc < s. Then it 
is easy to check that the equation X[z,fc]7r''= — Pi[i, fc]vr''= 
has exactly q''' solutions for X[i,k]. It follows that (|6]l has 
exactly gf"('l^ ^'■"s) solutions. Finally, by using the fact that 
J2k=i ^k = YlnZi i[i^i+i - i^t)y wc Complete the proof. ■ 

IV. Matrices under Row Constraints 

In this section, we study the class of matrices in P"^™ 
whose rows are constrained to be elements of P^. In particular, 
we will present two constructions and several counting results 
for this class of matrices. These constructions and results are 
of primary importance to our study of capacities and coding 
schemes in later sections. 

A. TT-adic Decomposition 

Let R^^i^'^ denote the set of matrices in P"><™ whose rows 
are elements of P^. Then the size of P"^*' is 



l^nxp, ^ |i?M|"^ n|A.| 



(7) 



since there are |P^| = q'f^' choices for each row. Taking the 
logarithm on both sides of ^, we obtain 



logJP"^^|=n|/i| 



(8) 



Every matrix X G P"^^ can be constructed based on its 
TT-adic decomposition 



X = Xo+ ttXi 



TI"^ Xs^l, 



with each auxiliary matrix Xi (i = 0, . . . , s — 1) satisfying: 

1) Xi[l:n, l:^i+i] is an arbitrary matrix over TZ{R, tt), and 

2) all other entries in Xi are zero. 

The construction is illustrated in Fig. [T] Clearly, this con- 
struction provides a one-to-one mapping from sequences of 
7i|/i| q-ary symbols to matrices in P"^''. 

B. Row Canonical Forms in 7^(P"x^) 

Let 7^(P"x^) denote the set of matrices in R"^!^'- whose 
shape is k. Then |7^(P"x^)| = unless n < n and n < jjl 
(written k < n, /i). The first constraint comes from the fact 
that the row canonical form of a matrix in P"^^' has at most n 
nonzero rows. The second constraint comes from the fact that 
row A is a submodule of P^, for any A G R"^^. Hence, we 
will assume that k ^ n,fi in the rest of this paper As we will 



see, the set 7^(_R"^^), together with the row canonical forms 
in 7^(_R"^^), plays a crucial role in our coding schemes. 

We now enumerate the row canonical forms in 7^(i?"^^). 
We need the following lemma. 

Lemma 1: There is a one-to-one correspondence between 
row canonical forms in 7^(_R"^^) and submodules of i?^ with 
shape K. 

Proof: Let S denote the set of row canonical forms in 
7^(_R"^^), and let Q denote the set of submodules of i?^ 
with shape k. Let (f) : S ^i- Q ht the map that takes a matrix 
B G 5 to its row module rowB. We will show that is a 
one-to-one correspondence. 

If 4){Bi) — 4'{B2) then Bi and B2 are left-equivalent, and 
so i?2 is a row canonical form of Bi and vice-versa. By the 
uniqueness of the row canonical form, we have Bi ^ B2', thus 
(j) is injective. 

Now let M be a submodule of i?^ with shape M = k, and 
construct a matrix A such that every element in M is a row 
of A. Clearly, row A = M and shape A = k. In particular, 
RCF(A) has Kg nonzero rows. Since k < n, we have Kg £ n. 
Let B be the submatrix of RCF(A) consisting of the top n 
rows. Then B contains all the nonzero rows of RCF(A), which 
means that shape i? = shape A — k. Hence, B e 7^(_R"^^), 
and the map is surjective. ■ 

By Lemma [T] the number of row canonical forms in 
7^(i?"^^) is [^] . It is helpful to bound this number as well 
as the logarithm of this number Combining (|5]l and the fact 
that 



q' 



k{m—k) 



< 



< 4^q 



k{m~k) 



(see, e.g.. 



Lemma 4]), we have 



-,12^=1 l^iit^i 



< 



< /^'iqI2t=l^^^i^^^-l^^ 



(9) 



Taking logarithms, we obtain 



i=l 



"^Kiifli - Ki) < log^ < ^Ki(^j - Ki) +Sl0gg4. 



(10) 



Example 11: Let R == Z4, and let n = 2, ^ = (2, 3), k = 
(1,2). Then by Lemma [T] there are 18 row canonical forms 
in 7^(i?"^^). These 18 row canonical forms can be classified 
into 4 categories based on the positions of their pivots; 
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The first category contains 8 row canonical forms, namely. 
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The second category contains 4 row canonical forms, namely. 
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I * * * * 
I * * * * 



]_ I ^ * * ^ 



-1^3- 



Fig. 2. Illustration of the construction of principal row canonical forms for 
TkCK"'^^) with s = 3, n = 6, A* = (4, 6, 8), and k = (2, 3, 4). 



The third category contains 4 row canonical forms, namely. 



ri oi 
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The fourth category contains 2 row canonical forms, namely. 



1 0' 




[2 1 0] 


2 




2 



Clearly, the first category contains a significant portion of all 
possible row canonical forms. ■ 

Motivated by the above example, we introduce principal 
row canonical forms that make up a significant portion of all 
possible row canonical forms in 7^(i?"^^). 

A row canonical form in 7^(i?"^^) is called principal if 
its diagonal entries di,d2, ■ ■ ■ ,dr (r = min{n, m}) have the 
following form: 

1 _s-l 



di 



,,7r 



,0,...,0. 



(11) 
Clearly, the first category in Example [TT] contains all principal 
row canonical forms for 7^(Z4^^) with n — 2, n — (2,3) 

and K= (1,2). 

Proposition 6: Every principal row canonical form X G 
Tk,{R"'^^) can be constructed based on its 7r-adic decomposi- 
tion 



X = Xo + ttXi + • 

with each auxiliary matrix Xi (i 
following conditions: 

1) Xi[l:Ki+i, l-.Ki+i] ^ diag(0, 



0, . . . , s — 1) satisfying the 
,0,1,...,1), 



2) Xi[l:Ki+i, Ki+i + l:/x,;+i] can be any matrix over 
n{R, vr), and 

3) all other entries in Xi are zero. 

This construction is illustrated in Fig. |2] Clearly, this 
construction provides a one-to-one mapping from sequences 
of X]i=i '^iif^i ~ ^i) 9-ary symbols to principal row canon- 
ical forms in 7^(i?"^^). Note that the number of principal 
row canonical forms in 7^(-R"^^) is q5i;i=i Ki(A'i-Ki)^ which 
is comparable to the number of row canonical forms in 
7;(i?"''^) in total. 

Proof: We will show that (i) every X constructed as above 
is a principal row canonical form, and (ii) every principal 
row canonical form has a 7r-adic decomposition following the 
above conditions. 

We begin with Claim (i). First, we track the diagonal entries 
in X. Clearly, by construction, the first ki diagonal entries in 
X are 1; they are contributed by Xq. The next K2 — K1 diagonal 



10 



entries in X are tt; they are contributed by Xi. Continuing 
this argument, we conclude that the diagonal entries in X are 
indeed of the form (fTTT i. 

Second, we show that X satisfies all the four conditions for 
row canonical forms. 

1) By construction, the first Kg rows of X are the only 
nonzero rows. Hence, X satisfies Condition 1. 

2) It suffices to show that the nonzero diagonal entries 
are precisely the pivots in X. Suppose that the ith 
diagonal entry X[i,i] = vr'. Then by construction, tt' 
is contributed by Xi and ki < i < Ki+i. Note that 
for each auxiliary matrix Xii, only the first kj'+i rows 
are nonzero. Thus, the jth row in Xii is zero for all 
I' = 0,...,/ - 1. In particular, Xi>[i,j] = 0, for all 
I' — 0, ... ,1 ~ 1 and for all j > i. Therefore, we have, 
for all j > i, 

s~l 

l'=Q 
s-1 

i'=i 

s-l 

i'=i 

That is, every X[i,j] is a multiple of tt' whenever 
j > i. On the other hand, by construction, X[i. j] = 
whenever j < i. It follows that X[i, i] is indeed the pivot 
of row i. Hence, X satisfies Condition 2. 

3) Since the nonzero diagonal entries are the pivots, X 
satisfies Condition 3. 

4) Suppose that the ith pivot X[z,i] = tt'. Then, we have 
Ki < i < K;+i. Note that for each auxiliary matrix Xi', 
all other entries in column i are zero as long as I' > I. 
Thus, we have, for all j ^ i, 



/-I 

= Y,7T'X,[J,^]. 



l'=0 

It follows that X[j,i] e 7^(i?, tt') for all j ^ i. Hence, 

X satisfies Condition 4. 
We turn now to Claim (ii). Let X be a principal row 
canonical form in 7^(i?"^'^). Then the diagonal entries in 
each Xi must satisfy 



Xi[l, 1], . . . , Xi[Ki+i, Ki+l] — 0, 



,0,1,, 



,1 



C. General Matrices in T^{W'^'') 

Here, we enumerate the matrices in 7^(_R"^^). We need 
the following lemma. 

Lemma 2: The number of matrices in R^^f^ having a given 
row canonical form in 7^(i?"^^) is equal to 



Iff' 



n (1 - q'-^- 



i=o 



Proof: Let _B be a row canonical form in 7^(i?"^^). Let 
B be the submatrix of B consisting of only nonzero rows. 
Clearly, B £ i?''=^''. We would like to count the number of 
matrices in R'^^f^ having the row canonical form B. 

By Proposition |5] every matrix A with RCF(A) — B has 
(7"^^=i ^(Ki+i-Ki) decompositions of the form A — CB for 
some full-column-rank C G i?"^"'. Hence, the number of 
matrices in R^^f^ having the row canonical form B is equal 
to the number of full-column-rank matrices of size n x k^ 
divided by q"^i=i '(Ki+i-^i)^ which can be simplified to 

1^ llli=0 {^-1 )■ ■ 

We can partition all the matrices in 7^(i?"^'') based on 
their row canonical forms: two matrices belong to the same 
class if and only if they have the same row canonical form. By 
Lemma [T] the number of such classes is [^] .By Lemma |2] 

the number of matrices in each class is \R"^'^\ nr=o ^^ ~ 
q*^"). Combining these two results gives us the following 
theorem. 



Theorem 3: The size of 7^(i?"^^) is given by 



\%iR'''"')\ = |i?"^'"| Yl (1 -?'"") 



?:=o 



(12) 



Theorem[3]counts the number of matrices in R"-^f^ of shape 
K. In particular, when the chain length s = 1, R becomes Fg, 
and this counting result becomes HtiO (l" ~1^) [k^] ' which 
is the number of n x ji^ matrices of rank ni. We note that 
Theorem |3] generalizes a theorem of [14] from square matrices 
to general matrices and from Galois rings to finite chain rings. 

Taking logarithms on both sides of (fT2] i. we have 



iogjr.(i?""^)l = iog. 



log,|i?""1 



log, WH-q'-). 



4=0 



Combining this with dS) and ( fTol i. we obtain 



Moreover, since X satisfies Condition 4, it follows that each 
Xi satisfies the first condition described above. Since X 
satisfies Condition 2, it follows that Xi[Ki+i + l:n, 1:to] is 
a zero matrix. Finally, due to the constraints imposed by /i, 
Xj[l:n, /Xi+i + l:m] is a zero matrix for all i. Therefore, each 
Xi satisfies the second and third conditions. This completes 
the proof. ■ 



Y^ n,{n + M* - i^i) + logg n (^ ~ «' ") 

1=1 i=0 

S Ks — 1 

^At,(n + ^, -«;,)+ log, Jl (l-g'-") + slog,4. (13) 

1=1 1=0 
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V. Matrix Channels over Finite Chain Rings 

Let i? be a {q, s) chain ring, and let fi be an s-shape. 
Consider a matrix channel given by 



Y = AX + BZ 



(14) 



where X G _R"^^ and F e R^^'^ are the input and 
output matrices, respectively. The error matrix Z £ i?*^^ 
and the transfer matiices A 6 ^JVxn ^^^^ ^ ^ i?^^* 
are random matrices with some joint distribution conditioned 
on X. Clearly, (fT4l i is an instance of the discrete memoryless 
channel (A',py|x,3^) with input alphabet X = i?"^'', output 
alphabet y = R^^'^ and channel transition probability py\x- 
The capacity of this channel is given by 

C = maxI{X;Y) 

Px 

where px is the input distribution. We assume that logarithms 
are taken to the base q, so the capacity is given in g-ary units 
per channel use. 

Note that we may think of the rows of X, Y and Z as 
packets over an ambient space i?^. To support this ambient 
space, the length of a packet, denoted by m, is equal to fig. In 
many situations, it is useful to understand the capacity scaling 
as the packet length grows. For that reason, we introduce a 
notion of asymptotic capacity 



C= lim 



1 



-C= lim 



1 



rC, 



7n~^co n\^\ 



m-foo n\IJ,\m 



where we assume that n = n/m and /2 — {fti , . . . ,p,s) = /i/m 
are fixed. Note that C is normalized such that C ~ I if the 
channel is noiseless (i.e., A = I and Z — 0). 

Let T be an s-shape such that t ^ t, fi. Similarly to lHJ, we 
are primarily interested in the case where the transfer matrix A 
is uniform over GL„ (i?) (in particular, N = n), B is uniform 
over 7i(i?"''*), Z is uniform over 7;(i?* ''''), and X, A, B and 
Z are statistically independent. (The model in fT| is recovered 
when R = ¥q, jj, = m, and r — t.) In this case, the channel 
model (O is equivalent to 



Y = A{X + A-^BZ) = A{X + W), 



(15) 



where A e GL„(i?) and W = A-'^BZ e Tr{R''''^) are 
chosen uniformly at random and independently from any other 
variables (as shown in Appendix |A]l. We will focus on the 
channel model ( fTSl l in the rest of this paper. 



VI. The Multiplicative Matrix Channel 

As a first special case, following LU, we consider the 
multiplicative matrix channel (MMC) defined by the law 

Y ^AX, 

where A is uniform over GL„(i?) and independent from X. 
This model is a special case of the channel model ( fTSl l with 
W ^0. 



A. Capacity 

The capacity of the MMC follows directly from Proposi- 
tion 1 of |1|, which we reproduce here for convenience. 

Proposition 7: fl\ Proposition 1] Let Q he a finite group 
that acts on a finite set S. Consider a channel with input 
variable X € S and output variable F G 5 given by y = AX, 
where A G G i& drawn uniformly at random and independently 
from X. The capacity of this channel, in bits per channel use, 
is given by 

c^iog^\s/g\, 

where \S/Q\ is the number of orbits of S under the action 
of Q. Any complete set of representatives of the orbits, each 
transmitted with equal probability, is a capacity-achieving 
code. 

Now note that the general linear group GL„(i?) acts on 
j^nxfj. ^y left-multiplication, and each orbit is the set of 
matrices having the same row canonical form. Hence, the 
number of orbits is precisely the number of row canonical 
forms in i?"^'^. By Lemma [1] the number of row canonical 
forms in R^'^t^ is equal to the number of sub modules of i?'' 
with shape ^ n, /i. Therefore, we have the following theorem: 

Theorem 4: The capacity of the MMC, in q-ary symbols 
per channel use, is given by 



C, 



MMC 



log? Y. 



A capacity-achieving code C C /^"xm consists of all possible 
row canonical forms in i?"^^. 

Theorem H suggests that the transmitter could encode infor- 
mation in the choice of submodules. This naturally generalizes 
the "transmission via subspaces" strategy in ||201 . 

Corollary 1: The capacity Cmmc is bounded by 



^Ki(^i-Kj) < Cmmc < ^ Ki(Mj-Kj) + log^4^ 

i=l j=l 



71 + S 



where Hi — min{n, [/ii/2j} for all i. 

Proof: First, since k = (ki, . . . , Kg) d: n, (x, we have 



(16) 



Cmmc = log^ ^ 

A ^ n , /^ 

> log, ^^' 



s 
> ^K.,(^i - Ki), 
i=l 

where the second inequaUty follows from (fTOl l. 
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Second, we have 

Cmmc = log. 






<log, E 4^ 

A^n./i 



^E.A.(^. 



<log« 



IDi Ki{lli-Ki 



<log,4 



n - 



lEiKilA'i-Ki) 



^Kj(/ii - K,,) +l0gg4' 



i=l 

where the first inequaUty follows from (fTOt . the second in- 
equality follows from the fact that k maximizes the quantity 
Tlii ^i{f^i~^i) subject to the constraint \<n,^, and the third 
inequality follows from the fact that the number of shapes 
satisfying A ^ n,/i is upper-bounded by ("^'')- ■ 

We next turn to the asymptotic capacity of the MMC. 

Theorem 5: The asymptotic capacity Cmmc is given by 

n\n\ 

where R = K/m with Ki — min{n, [/ii/2j} for all i. 

Proof: This follows from Corollary [T] and the fact that 
^ log„ 4" ("+") ^ 0, as TO ^ cx). ■ 

Theorem |5] implies that the shape k given by Ki — 
minjri, [/ii/2j} (1 < J < s) is "typical" among the shapes of 
all possible row canonical forms in i?"^**. In other words, the 
row canonical forms of shape k make up a significant portion 
of all possible row canonical forms. Hence, the transmitter 
may encode information in the choice of row canonical forms 
of shape k instead of all row canonical forms. 

B. A Simple Coding Scheme 

In this section, we present a simple coding scheme that 
achieves the asymptotic capacity in Theorem |5] The key 
idea is to make the codebook the set of all principal row 
canonical forms for 7^(_R"^^). In other words, we employ 
two "reductions" in the code construction. First, we move from 
all row canonical forms in R'^'^t^ to all row canonical forms 
in 7^(i?"^^), as suggested by Theorem |5] Then, we move 
from all row canonical forms in 7^(_R"^^) to all principal row 
canonical forms in 7^(_R"^^). With these two reductions, our 
coding scheme not only achieves the asymptotic capacity, but 
also admits fast encoding and decoding. 

1) Encoding: The input matrix X is chosen from the set 
of principal row canonical forms for Tk{R^^^) by using the 
construction presented in Section IIV-BI Clearly, the encoding 

rate of the scheme is i?MMC = X]i=i ^i(A*i ~ i^i)- 

2) Decoding: Upon receiving Y = AX, the decoder simply 
computes the row canonical form of Y. The decoding is always 
correct by the uniqueness of the row canonical form. By 
comparing the encoding rate with the asymptotic capacity, we 
have the following theorem. 

Theorem 6: The coding scheme described above achieves 
the asymptotic capacity (fTTI l. 



VII. The Additive Matrix Channel 

In this section, we consider the additive matrix channel 
(AMC) defined by the law 

Y = X + W, 

where W is uniform over Tt{R"'^'^) and independent from 
X. This model is a special case of the channel model (flST l 
with A^ I. 

A. Capacity 

Theorem 7: The capacity of the AMC, in g-ary symbols 
per channel use, is given by 

Camc = log, |i?""^| - log, |r.(i?"^^)|, 

achieved by the uniform input distribution. 

Proof: The AMC is an example of a symmetric discrete 
memoryless channel, whose capacity is achieved by the uni- 
form input distribution. Note that when X is uniform over 
_R"^^, so is Y. Thus, we have 

Camc - H{Y) ~ H{W) = log, |i?"x^| - log, |r.(i?"^^)|. 

■ 

Corollary 2: The capacity Camc is bounded by 

S Ts-l 

^(n - n){n, n) - log, 4« n (1 - '?'"") 

i=l i=0 

< Camc < 

S Ts — 1 

J2{n - nKi^r - n) - log, n (1 - '?'"")• (18) 



i=0 



(O. 



Proof: It follows immediately from Theorem Q and (|8]l. 



We next turn to the asymptotic behavior of the AMC. 
Theorem 8: The asymptotic capacity Camc is given by 

Y.Uli^ - ^i)il^i -^i) 



Camc — 



n|^| 



(19) 



Proof: It follows from Corollary |2] and the fact that 



^iog,4^n(i-,-)^o, 



as ?77, — >■ oo. 



B. Coding Scheme 

We focus on a special case when t — t, and present a 
coding scheme based on the idea of error-trapping in [1 1. This 
scheme achieves the asymptotic capacity for this special case. 

1) Encoding: Set v > t. The input matrix X is constructed 

as , , 



U ' 



X ^ 



where the size of U is {n — v) x (to, — w), and the sizes of other 
zero matrices are chosen to make X an nxm matrix. Here, U 
is chosen from the set Ri'^^^)^it^^^) by using the construction 
in Section |IV] (as illustrated in Fig. [3]l. Clearly, the encoding 
rate of the scheme is Ramc = X]i=i("- ~ ^)(Pi ~ ^)- 
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^0 = 



^IH^ 



Xi = 









-M2- 



X.,= 






Fig. 3. Illustration of the AMC encoding scheme for 
At = (4, 6, 8), and v = 2. 



-1^3- 



3, n 



2) Decoding: Following [l], we write the noise matrix W 



W = BZ 



Bi 

B2 



[Zi ^2], 



where Bi e i?"^*, B2 G i?(»-^')xt^ Zi e R^""" and Z2 £ 
■pjtx(m-v)^ The received matrix Y is then given by 



Y = X + W = 



BiZi 
B2Z1 



B1Z2 
U + B2Z2 



Similar to HI, we define that the error trapping is successful 
if shape i?iZi — t. Assume that this is the case. Then by 
Proposition [T] 3, we have shape i?i = shape Zi — t. Consider 
the submatrix consisting of the first v columns of Y. Since 
shape BiZi ~ t, the rows of B2Z1 are completely spanned 
by the rows of BiZi. That is, rowi32Zi ^ rowi?iZi. Thus, 
there exists some matrix T such that B2Z1 = TBiZi. Since 
Zi is full row rank, by Theorem |2] B2Z1 = TBiZi implies 
B2 = TBi. It follows that 



T 



Bi 
B2 



Bi 





where T 



Note also that TX = X. Thus, 



TY ^TX + TW 



BiZi 




/ 
-T I 



B1Z2 
U 



from which the data matrix U is readily obtained. 

The decoding is summarized as follows. The decoder ob- 
serves BiZi, B1Z2, and B2Z1 thanks to the error traps. The 
decoder then checks the condition shape BiZi ~ t. If the con- 
dition does not hold, the decoder declares a failure. Otherwise, 
the decoder finds a matrix T such that B2Z1 = TBiZi (which 
means B2 — TBi). Since B2 = TBi, the decoder can recover 
B2Z2 by using the relation B2Z2 — TBiZ2- Clearly, the error 
probability of the scheme is zero. The failure probability of 
the scheme is 

Pf =Pr [shape Bi^i ^ t]. 

Lemma 3: The failure probability Pf of the above scheme 
is upper-bounded by Pf < t^t^t ■ 



Proof: If Bi and Zi are full rank, then shape i?iZi = t. 
Hence, by the union bound, the failure probability 

Pf < Pt[Zi is not full rank] + Pr[Bi is not full rank]. 

Now consider the probability that Zi is full rank. Recall that 
Z ^ R^^'^ is a full -rank matrix chosen uniformly at random. 
An equivalent way of generating Z is to first generate the 



entries of a matrix E G _R*^'' uniformly at random, and then 
discard E if it is not full rank. This suggests that 

Pi-[Zi is full rank] = Pr[£'i is full rank | E is full rank] 
> Pr[£'i is full rank], 

where Ei consists of the first v columns of E. Thus, 



Pr[Zi is full rank] > |7;(i?*''")|/|i?*''' 
t-i 

i=0 

t-1 

= 11(1-9'"") 

4=0 

> 1 



)/q' 



t 



q' 



Similarly, we can show that 

Pr[Bi is full rank] > 1 



Therefore, the failure probability Pf < rj*_t . ■ 

Recall that the encoding rate of the scheme is i?AMC = 
jy^-iin — v){fi.i — v). Thus, if we set v such that 

v-t 



v-t 



cx), and 



as m ^> 00, then we have Pf -> and Ramc = -^^f^ -^ 
Camc- Therefore, we have the following theorem. 

Theorem 9: The coding scheme described above can 
achieve the capacity expression (fT9] l for the special case when 

T = t. 

VIII. The Additive-Multiplicative Matrix Channel 

In this section, we consider the additive-multiplicative ma- 
trix channel (AMMC) defined by the law 

Y ^A{X + W), 

where A G GL„(i?) and W G Tr{R'^'''') are uniformly 
distributed and independent from any other variables. 

A. Capacity Bounds 

Theorem 10: The capacity of the AMMC, in q-sry symbols 
per channel use, is upper-bounded by 

CamMC < CmmC - log, |r.(i?"^^)| + log, Y. ITr'l-R"''")!- 

(20) 



t' -<T 



Proof: First, using the same information-theoretic argu- 
ment as in |1), we have 



I{X;Y)<CMMC~log„\Tr{R'' 



XpN 



H{W\X,Y) 



Next, we upper bound the term H{W\X,Y). Let k — 
shape y. Let S be the Smith normal form of Y. Then 
S contains Kg nonzero diagonal entries. Thus, Y can be 
expressed as 



Y = PSQ = [Fi P2] 



Sn 




Qi 
Q2 



— PiSuQi, 
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where Pi e 7?"^"^% Qi e i?''=x'", and 5ii e P'^^x"-. 
Note that 

X + VK = A^^r = A*Qi, 

where A* ~ A^^PiSu. Since Qi consists of the first k^ rows 
of an invertible matrix Q, Qi is a full-rank matrix. In particu- 
lar, Qi contains an invertible Kg x Kg submatrix. By reordering 
columns if necessary, we can assume that the left Kg x Kg 
submatrix of Qi is invertible. Write Qi = [Qn Q12], 
X = [Xi X2] and W = [Wi W2], where Qn, Xi, and 
Wi have Kg columns. We have 

A* = (Xi + Wi)Q^,^ and W^2 - A*Qi2 - X2. 

This suggests that W2 can be computed from Wi if X and Y 
are known. Thus, 

i7(VF|X,r) ^H{Wi\X,Y) <i7(T^i|shapey). 

Since Wi is an n x Kg matrix with shape Wi < r, we have 

H{Wi\shapeY = k) < log^ ^ |r.,(i?"^''OI, 

which is maximized when k,s — n. Hence, 

i/(W^i|shaper) < log^ Y. \Tr'{R"'''')\. 

It follows that H{W\X,Y) < log^ Er'^r l'^'(-R"''")l' 
which completes the proof. ■ 

Corollary 3: The capacity Cammc is upper-bounded by 

s s 

Cammc < ^if^i - Ci)6 + X!^" ~ '^'^'^^ + Sslog^ 4 
1=1 i=i 

+ log, ("r) + log, (^r) - log, n (1 - 9^'")' 

i=0 

where ^^ — min{n, [/ii/2j} for all i. In particular, when /i >^ 
2n, the upper bound reduces to 

s 

Cammc < ^(»^ - n){m -~n) + 2s log^ 4 

i=l 

+ log, (T) + log, i^'t') - log, n\l - q'-n- 
Proof: By ( fTSI l. we have 

CmMC <Yi^J■^- ^t)^t + S log^ 4 + log^ f 

1=1 ^ ^ 

By ( fT3] l, we have 

s 

^ log, ir.(P""^)| < - ^(n + M,; - n)n 

Ts-1 

-log, Ha -9"") 



i=0 



Note that 



< 4SQELl(2n— rOn' 



where the first inequality comes from (fTZt . and the second 
inequality comes from ^ and ( fTOl i. Hence, 



-j-Z^-p 



t' ^T 



where the second inequality comes from the fact that r 
maximizes the quantity q^l=\^^'^^'^i)'^i and the fact that the 
number of shapes r' with r' ^ r is upper-bounded by i^'^^') ■ 
Therefore, we have 

log, Y. ir.'(i?""")l<E(2"-^^)^*+^i°g.4+iog,rr). 

t' ^T ?'=1 

Combining all the above results, we have obtained the upper 
bound. In particular, when /i ^ 2n, we have ^,; = n for all i. 
Substituting this into the upper bound completes the proof. ■ 
We next study the asymptotic behavior of Cammc- 

Theorem 11: When /i ^ 2n, the asymptotic capacity 
Cammc is upper-bounded by 



Cammc < 



E»=i("~'r»)(A^» 
n\[i\ 



(21) 



Proof: This follows from Corollary |3] 



B. A Coding Scheme 

We again focus on the special case when t ^ t. We 
describe a coding scheme that achieves the asymptotic bound 
in Theorem [TT] when n ^ 2n for this special case. 

1) Encoding: The encoding is a direct combination of the 
encoding strategies for the MMC and the AMC. Specifically, 
with V > t, the input matrix X is constructed as 



X 



q 

X 



where the size of X is {n — u) x (m — v), and the sizes 
of other zero matrices are readily available. Here, X is 
chosen from the set of principal row canonical forms for 
j-^(^jl{n-v)x{f,-v)^ by using the construction in Section H^ll 
where Ki = min{n— w, [(/ii — u)/2j } for all i. The encoding is 
illustrated in Fig. |4] Clearly, the encoding rate of the scheme is 
Pammc = Y^i=i i^iit^i -V ~ Ki). In paiticular, when /i >z 2n, 
we have [(/i.^ — v)/2\ > n — v for all i. Thus, Ki ^ n — v for 
all i, and the encoding rate is Pammc = X]i=i('^^^)(/"i ^")- 






Xu- 









n*] 



X, 






X, 



* * * 

* * * 
1 * * * 



^)H^ 



-M2- 



-fki- 



Fig. 4. Illustration of the AMMC encoding scheme for s = 3, n 
1, = 2, ^ = (4, 6, 8), so that k = (1, 2, 3). 
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2) Decoding: The decoder receives Y = A{X + W) and 
attempts to recover X from the row canonical form of Y. 
Again, we write the noise matrix W as 



W = BZ = 



Bi 
B2 



From Section IVIII we have 

X + W = 



BiZi 
BiZ\ 



[Zi Z2] 



_ B\Zi 
X + B2Z2 



Following [1|, we define error trapping to be successful 
if shape BiZi = t. Assume that this is the case. From 
Section IVIII there exists some matrix T G GL„ (i?) such that 



BiZi 




BIZ2 




\Bi 


0' 




\Zi 


Z2] 


X 







/ 







X 



T{X + W) = 
Note further that 

RCF 



for some Zi € i?*^^ in row canonical form and some Z2 G 
j^tx{m-v)_ It follows that 



\Zi 


Z2] 


\ 


\Zi 


Z2 





X 


)' 





X 



RCF(X + W)^ RCF 



^r^TT 


/ 


B 


1 0] 




\Zi 


Z2] 


r^i 


V 

2 


C 

7 

'2 


I 







X 





X 











D_ 











where the bottom v — t rows are all-zeros. 

Since A is invertible, we have RCF(y) = RCF(X + W), 
from which X can be readily obtained. Hence, decoding 
amounts to computing the row canonical form, whose com- 
plexity is 0{nmT[vm{n,ra}) basic operations over R. 

The decoding can be summarized as follows. First, the 
decoder computes RCF(F). Second, the decoder checks the 
condition shape BiZi = t. If the condition does not hold, the 
decoder declares a failure. Otherwise, the decoder outputs X 
from RCF(r). 

Let Y denote the left-most n columns of RCF(y), i.e., 
Y = RCF(y)[l:ri, V.n]. We note that shapcBiZi = t if and 
only if shape Y — t + n. Hence, the error probability of the 
scheme is zero, and the failure probability Pf of the scheme 
is bounded by Pf < i^t-t (as shown in Section IVIlb . 

Recall that the encoding rate i?AMMC — X]i=i('^^''^)(/^i~") 
when /i >z 2n. Thus, if we set v such that v — t ^ 00 and 
^^^ — > 0, as TO — >■ cx), we have Pf -^ 0, and /Jammc = 
4t ^ ^'-^<'lT*'](^'-''^) for fi >- 2n. 

Theorem 12: When t — t and /i >r 2n, the coding scheme 
described above can achieve the upper bound (l2Tt . 

Finally, we would like to comment on the condition of r = t 
and /i >r 2n. Let Zi be the submatrix of Z consisting of the 
leftmost fii columns. We notice that the condition r = t is 
equivalent to the condition Zi G _R*^^i is full rank. According 
to Lemma [3] when the error matrix Z is uniform over i?*^'^, 
the probability that Zi is full rank is lower bounded by 

Pr[Zi is full rank] > 1 



Clearly, when /ii is sufficiently large, this probability tends to 
1. We also notice that the condition fi ^ 2n is equivalent to 
the condition /ii > 2n, which holds when /ii is large enough. 
To sum up, the condition of r = t and fi ^ 2n seems to be 
mild in the scenario where /ii can be made sufficiently large. 
Such a scenario can be found in some existing constructions of 
nested-lattice-based physical-layer network coding (see, e.g., 
im Examples 8 and 9]). Nevertheless, the general scenario stiU 
remains open. 

IX. Conclusions 

In this work, we have studied the matrix channel Y = 
AX + BZ over a finite chain ring. Under the assumption 
that A is uniform over all invertible matrices and BZ is 
uniform over all matrices of shape t, we have derived tight 
capacity results and provided polynomial-complexity capacity- 
achieving coding schemes, which naturally extend the work 
of jTj from finite fields to finite chain rings. Our extension 
requires the use of row canonical forms, as well as several new 
enumeration results and construction methods, for matrices 
over finite chain rings, which may be of independent interest. 

We believe that there is still much work to be done in 
this area. One direction would be to relax the assumption 
on A and BZ. Following this direction, we have explored 
a particular case when A can be any matrix and BZ = in 
II2TI . Another direction would be to find other applications of 
the row canonical form developed in this paper 

Appendix 

A. Proof of the fact that W = A^^BZ is uniform over 
7^(_R"^^) and independent of A. 

First, we will show that W is uniform over Tt{R"'^^)- We 
need the following lemma. 

Lemma 4: The map ip : 7;(i?"^*) x 7;(i?*^^) -^ 
7;(i?"^^) defined by 

ifiA, B) = AB, VA e TtiR"'"'), B e Tr{R''"') 

is surjective. Furthermore, the pre-image (^^^(C) of every 
matrix C G %-{R^^^) has the same size. 

Proof: First, we will show that the map tp~^ is surjective. 
For any C G Tt{R^^^), let C be the row canonical form of 
C. Then C ^ PC for some invertible P G GL„(i?). Let 
C be the submatrix of C consisting of the top t rows. Let 
P = [Pi P2], where Fi G i?"""* and P2 G i?"x ("-*). Then 
we have 



C = PC= [Pi P2] 



= PiC. 



q 



l + A"! 



Since Pi is full-column-rank and C has shape r. Pi G 
TtiR"''*) and C G TriR*'"'')- Hence, cp is surjective. 

Next, we will show that \(f~^{Ci)\ ~ \(p^^{C2)\ for any 
two matrices Ci,C2 G Tt{R^^^)- Since Ci is equivalent to 
C2, there exist invertible matrices P and Q such that C2 — 
PCiQ. Let {Ai,Bi) be a pre-image of Ci, i.e., AiBi = 
Ci. We will show that {PAi,BiQ) is a pre-image of C2. 
Since PAiBiQ — PCiQ — C2, we only need to show that 
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PAi e TKi?"""*) and BiQ e TriR^'"''). Since P preserves 
the row span, we have PAi e Tt{R"^'^)- Since PAi is full- 
column-rank, we have tow BiQ — rowC2. It follows that 
BiQ e Tt{R*'^^')- Hence, {PAi,BiQ) is indeed a pre-image 
of C2. We have |(p~^(Ci)| < \tp-^{C2)\. Similarly, we can 
show that |iy9^^(C2)| < |iy9^^(Ci)|. This completes the proof. 

■ 

By Lemma 21 if i? is uniform over 7t(-R"^*) and Z is 
uniform over Tt{R*^^), BZ is uniform over Tt{R"'^^)- 
Moreover, if A^^ is uniform over Tn{R"^^) and BZ is 
uniform over 7;(i?"''^), A-^BZ is uniform over TriR"'^^'). 
That is, M^ is indeed uniform over Tt{R"'^^)- 

Second, we will show that W is independent of A. Note 
that for any given A, A^^BZ is uniform over Tt{R^^^^)- It 
follows that W is independent of A. 
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